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Chaos assisted instanton tunneling in one dimensional perturbed periodic potential 
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For the system with one-dimensional spatially periodic potential we demonstrate that small peri- 
odic in time perturbation results in appearance of chaotic instanton solutions. We estimate param- 
eter of local instability, width of stochastic layer and correlator for perturbed instanton solutions. 
Application of the instanton technique enables to calculate the amplitude of the tunneling, the form 
of the spectrum and the lower bound for width of the ground quasi-energy zone. 
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INTRODUCTION 

Tunneling as inherently quantum phenomenon attracts 
much attention . Its connection with classical chaos in 
semiclassical regime has also been discussed |2j, |3|- A 
number of works were devoted to semiclassical chaos as- 
sisted tunneling between symmetry related KAM-tori in 
systems with mixed dynamics (well developed chaotic re- 
gion coexists in phase space with regular islands) HUES- 
To describe chaos assisted tunneling in systems with 
mixed dynamics multi-level model Hamiltonian, primar- 
ily proposed in 0, is often used Q. Less attention has 
been payed to semiclassical tunneling in KAM systems 
(chaotic region is not widespread) Another way to 
describe semiclassical tunneling is based on solutions of 
Hamilton equations in imaginary time and path integral 
formalism M. Instanton technique [jj was used in a very 
few works 

In this work we consider one-dimensional quantum sys- 
tem with periodic in space potential affected by small 
periodic in time perturbation. We use methods created 
to describe chaos in classical Hamiltonian systems to in- 
vestigate essentially quantum phenomenon of tunneling. 
It is achieved in the framework of instanton technique, 
where solutions of Euclidian equations of motion (instan- 
tons) play dominating role, by the use of standard meth- 
ods from the viewpoint of chaos For the systems 
with periodic in time perturbation energy is no more 
an exact integral of motion and the language of quasi- 
energies is more adequate 1.2j. For some estimations 
energy as an adiabatic invariant can also be used [lflj . 
We study properties of chaotic instanton solutions and 
calculate the form of the spectrum and the lower bound 
for the width of the ground quasi-energy zone. 

Hamiltonian of the system under consideration is taken 
in the form 



H = -p 2 + lu 2 cos x - 
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unit. Considered cosine potential corresponds nonlinear 
oscillator. Phase space of nonlinear oscillator has the 
topology of cylinder (points (q,p) and (q + 2ir,p) are 
identified). Thus its (quasi-) energy spectrum is discrete. 
Chaos assisted tunneling between two major resonance 
islands (inside single potential well) for driven nonlinear 
oscillator has been studied numerically in |4(. In this 
work we study the system with x varying from — oo to 
+oo. This results in band structure of the (quasi-)energy 
spectrum 0] . Systems with spatially periodic potential 
were studied in the instanton physics |15j. Perturbation 
used in Q was exp loited in the systems exhibiting quan- 
tum chaos 0, 

There are papers devoted chaos assisted tunneling 
where some analytical predictions for billiard systems 
basing on the path integral formulation of quantum me- 
chanics are made Distinguishing feature of our work 
is analytical predictions for the system with smooth po- 
tential exploring for this purpose instanton technique 
adopted from quantum field theory 





ANALYSIS OF CHAOTIC INSTANTON 
SOLUTIONS 

For applying instanton technique we consider solutions 
of classical equations of motion in imaginary (Euclidian) 
time. Hamiltonian has the same form (translated on 
7r) in Euclidian time as in real one. 

Euclidian Hamiltonian of the system is H = Ho + V , 
where 



and 



H = \^p 2 - CJgCOSX, 



+ 00 

V = otTx S(r — nT). 



(2) 



(3) 



T is the real time period of perturbation, e describes the 
strength of perturbation. The mass of the particle equals 



Here Hq is nonperturbed Euclidian Hamiltonian of the 
system and V is the Euclidian potential of the pertur- 
bation. We also introduced coupling constant a <C 1 
instead of e = aT in order to simplify formulas. 
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Nonperturbed instanton solution describes the motion 
on the separatrix of the Hamiltonian (J2J) . It is well known 
that this separatrix is destroyed by any periodic pertur- 
bation and on its place stochastic layer is appeared [l3| . 
Perturbed instanton solutions correspond to the motion 
in vicinity of the separatrix inside the layer. Therefore in- 
stead of one instanton solution connecting neighbor max- 
ima of nonperturbed Euclidian potential (classical vac- 
uum states in real time potential) we obtain a manyfold 
of instanton solutions of Euclidian equations placed in- 
side the stochastic layer. 

We calculate parameter of local instability, width of 
stochastic layer and correlator for perturbed instanton 
solutions. It is convenient to describe dynamics of the 
system in action-angle variables . Equation of motion 
for action variable has the form 



+ 30 



I = — 



2 cos (rnvr) 
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Here u>(I) = dH^/dl is the nonlinear frequency 01 . In- 
stead of angle variable we introduce phase of external 
force ip defined by the relation ip = v = 2ir/T [l3|. 
Let H s = ujq to denote the energy of nonperturbed 
system on the separatrix. Continuous equations of mo- 
tion for / and ij) can be reduced to discrete mapping for 
the phase of external force in the vicinity of separatrix 

i\H-H s \<i) msini 



ip n+ i = ip n + B„ + K sinV>™, 



(5) 



where 



K 



uj \H-H S 



B n are some functions of H which exact form is not es- 
sential for our purposes. We assume following [13] that 
due to small value of perturbation the energy practically 
does not change with time and equals the energy of the 
nonperturbed system. The map 10 with arbitrary pa- 
rameter Kq was studied by many authors, for instance 
|l(j|. Particulary, it is known that at Kq > 1 motion 
is locally unstable and chaotic, whereas at Kq < 1 it is 
stable and regular. Thus Kq is the parameter of local 
instability. Condition Kq ~ 1 enables us to calculate the 
width of stochastic layer 



\H.-H b \ = 



CJ 



(6) 



here H{, is the energy value on the bound of stochastic 
layer. 

To calculate correlator for perturbed instanton solu- 
tions we use standard technique 18]. For the map Q 
correlator is 

i r 2 " 

R{t,t ) =— / d0 o exp{i(V>(r) - Vo)} 



exp 



T -Tp 
TR 



(7) 



here ipo = - 0( r o) and the time of correlations decay is 
tr = 2ir/(uj\nKo). Exponential decrease of correlator 
shows that dynamics of the instanton solutions inside the 
stochastic layer {Kq > 1) possesses the property of mix- 
ing (chaos) yjj . 

Note that perturbed one-instanton solution due to 
stochastic layer connects not only neighbor vacua of real 
time potential but also two arbitrary chosen vacua. Re- 
mind that to describe tunneling between non-neighbor 
vacua of nonperturbed system one has to take into ac- 
count contribution of multi-instanton configurations |l5| . 



THE CALCULATION OF THE TUNNELING 
AMPLITUDE AND GROUND ZONE WIDTH 

Let us consider the tunneling between neighbor vacua 
(from x Rj — 7r to x ~ 7r for distinctness) in presence of 
perturbation (jSJ)- In Euclidian time this tunneling pro- 
cess for the nonperturbed system is described by the so- 
lution of Euclidian equations of motion with asymptotes 
x = — 7r, p = at r = — oo and x = 7r, p = at 
t = +oo. There is only one solution satisfying these 
conditions for nonperturbed system © (one-instanton 
solution) 



V-To) 



4arctane" o(T " 



(8) 



Its Euclidian action is S mst — 8^o- The instanton's 
position is denoted by tq. Due to Euclidian equations 
of motion and anti-symmetry of x nst when time is in- 
verted with respect to the point To perturbation does 
not change the Euclidian action of the one-instanton so- 
lution |JSJ in the first order on the coupling constant 
S™* = S inst +0{a 2 ). The only manifestation of the 
perturbation in this approximation is the appearance of 
a number of the additional solutions of Euclidian equa- 
tions of motion with energies close to the energy of non- 
perturbed one-instanton solution and placed inside the 
stochastic layer. 

Let us consider firstly nonperturbed system at arbi- 
trary energy — ujq + e, < e < 2ujq. One half of trun- 
cated instanton action can be easily calculated 



o(e) 



S[x inst (T,e)} = 



2(o7q cos x — (— to 2 + e))dx 



-a(e) 



AJAco 2 -2eE o(e), ■ 



where ±a(e) = ±arcsin,/l 



2^ 



(9) 



are turning points, 



function E is the elliptic integral of the second kind. 
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Then tunneling amplitude in perturbed system can be 
found by integration over energy of the tunneling ampli- 
tude in nonperturbed system with the action ^ 

AH x(t)ritt 

A= J de J Dxex-p[-S[x mst {T,e)}] , (10) 

x(t)& — 7T 

where AH = 2\H S — Hb\ is the stochastic layer width. 
The contribution of the chaotic instanton solutions is 
taken into account by means of integration over e. Ex- 
pression l|lUfl shows that the probability of tunneling 
(square of the absolute value of the tunneling amplitude) 
grows while chaotic region spreads (AH increases). 

The result is obtained in the first order on coupling 
constant a and does not take into account possible struc- 



ture of stochastic layer is valid if the layer is narrow 
and is in agreement with results of numerical [2(| and 
real [f| experiments for similar problems. We have also 
correspondence in (|lf)|> with the nonperturbed case |l5l |. 
Namely, if a = then AH = and the single solution 
describing the motion on the separatrix (nonperturbed 
one-instanton solution) contributes to the tunneling am- 
plitude. 

Formula (|10|l can be made more transparent if we use 
the approximate form of action © at e < 2ujq 

S[x mst (T,e)}^8ujo-—. (11) 

LUq 

Then in Gauss approximation we obtain the following 
expression for the tunneling amplitude 



AH 



A 



de / dc oy /S[ai inst (T,e)]ex.p(-S[x in3t (T,e)]) a e - si ""Vs inat rF 



-Suio, 



(12) 



where integration over cq gives the contribution of zero 
modes, T is a time of the tunneling. Formula (|12|) up to 
the factor i has the same form in real Minkovski time. 
Expression (|12|) can be interpreted in the following 

way. Factor F = exp (^^p) > 1 in O differs per- 
turbed and nonperturbed amplitudes and includes the 
contribution of the layer. One can think about F as a 
number of instanton solutions inside the stochastic layer. 

We can find the minimal number of instanton solu- 
tions inside the stochastic layer. One can not distinguish 
instanton solutions within energy interval AE ~ 1/Ar 
(Heisenberg uncertainty relation). Here At 



de- 
notes the time interval of observation. Thus the energy 
interval between neighbor instantons is Ae ~ luq. There- 
fore parameter F can be found as follows 



F ~ 1 



AH 
~As~ 



(13) 



where unit takes into account the nonperturbed separa- 
trix. 

The form of the spectrum of the lower quasi-energy 
zone is obtained using the amplitude (|12fl by means 
of standard technique (multi-instanton contributions are 
taken into account) |l5j: 

Ee w Iruo - 2e~ s ' nst ^S mst F cos 6, (14) 



here continuous variable 9 parametrizes levels of the 
ground quasi-energy zone. Zone width 



AE « Ae- s '" st VS mst F 



(15) 



differs from the case of nonperturbed case by factor F 
reflecting the influence of perturbation. 



CONCLUSION 

We applied theory of classical chaos for investigation of 
the chaos assisted tunneling in terms of path integral for- 
malism in imaginary time and instanton technique. We 
found the parameter of local instability and the width of 
the stochastic layer. Exponential decrease of the correla- 
tor for any perturbed instanton solution was also demon- 
strated, that means it to possess the property of mix- 
ing. Then properties of the stochastic layer and classi- 
cal chaotic solutions in Euclidean space (chaotic instan- 
tons) were used for the calculation tunneling amplitude, 
ground quasi-energy zone spectrum in the presence of the 
perturbation and the zone width. 

General tendency for chaos assisted tunneling regime 
(in average - if we abstract from fluctuations) is the 
increase of tunneling amplitude (probability) as the 
strength of perturbation increases jy, |2fJ . It is confirmed 
here. The reason is the growth of the width of the chaotic 
layer and therefore the increase of the number of paths for 
particle to travel from one regular region to another. For 
small energies in Gauss approximation tunneling ampli- 
tude is increased by the factor F > 1. The life-time of the 
particle in the certain vacuum of the system decreases. 
It is connected with the widening of the (quasi-)energy 
zone <|15[) . 
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We would like to emphasize that obtained results are 
not consequences of the particular choice of the non- 
perturbed Hamiltonian J2J or perturbation J3J). Qual- 
itatively they are valid for more general class of one- 
dimensional nonpcrturbed Hamiltonians with quadratic 
dependence on momentum and spatially periodic poten- 
tial with single well in each period, as well as time de- 
pendence of homogeneous perturbation can be realized 
by any time periodic function. The reason is the univer- 
sality of the separatrix destruction mechanism in these 
potentials affected by time-periodic perturbation |l3| . 

Tunneling plays an important role in gauge field the- 
ories (instanton physics ||). Experimental discovery 
of QCD instantons for example is an important prob- 
lem [2lJ. Moreover it is known that classical gauge field 
theories are inherently chaotic Therefore the study 
of chaos assisted instanton tunneling in gauge field theo- 
ries based on chaos criterion in quantum field theory [2^] 
is also of essential interest. 
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